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Abstract 

Quadratic permutation polynomial interleavers over integer rings have recently received attention 
in practical turbo coding systems from deep space applications to mobile communications. In this 
correspondence, a necessary and sufficient condition that determines the least degree inverse of a 
quadratic permutation polynomial is proven. Moreover, an algorithm is provided to explicitly compute 
the inverse polynomials. 

Index Terms 



Turbo code, interleaver, algebraic, permutation polynomial, quadratic polynomial, inverse polyno- 
mial, zero polynomial, null polynomial. 



JAN. 23, 2011 



2 



I. Introduction 

Interleavers for turbo codes have been extensively investigated flTJ — [J3]| . Today the focus on 
interleaver constructions is not only for good error correction performance of the corresponding 
turbo codes but also for their hardware efficiency with respect to power consumption and speed. 
The work in [|2] opened the door to a class of polynomial based interleavers. In particular, 
quadratic permutation polynomials (QPP) were emphasized because of their simple construction 
and analysis. Their performance was shown to be excellent El, 0. The practical suitability of 
QPP interleavers has been considered in a deep space application [7] and in 3GPP long term 
evolution (LTE) lfl9ll . 

The inverse function for a QPP is also a permutation polynomial (PP) but is not necessarily 
a QPP 0. However, there exists a simple criterion for a QPP to admit a QPP inverse 0. A 
simple rule for finding good QPPs has been suggested in [3J. Some examples in do not have 
QPP inverses. Most of QPP interleavers proposed in 3GPP LTE [[T9l admit a quadratic inverse 
with the exception of 35 of them. 

In flU, a necessary and sufficient condition that determines the least degree inverse of a QPP by 
using Chinese remainder theorem and presenting the inverse function as a power series is given. 
As an example, an exact formula that determines the degree of the inverse PP is shown when 
the degree is no larger than 3. 

In this correspondence, we provide a necessary and sufficient condition by using linear congru- 
ence approach in [0 pp. 24-40] that determines the degree of the inverse when the degree is 
no larger than 50. The condition is characterized by an exact formula and consists of simple 
arithmetic comparisons. We further provide an algorithm to explicitly find the inverse PP(s). The 
algorithm is suitable for implementation since it consists of solving linear congruences. 
This correspondence is organized as follows. In section II, we briefly review PPs flTOll — [fT4ll over 
the integer ring Zjy and relevant results. The main result is derived in section III, and examples 
are given in section IV. Finally, conclusions are discussed in section V. 

II. Permutation Polynomial over Integer Rings 
In this section, we revisit the relevant facts about PPs and other additional results in number the- 
ory to make this paper self-contained. Given an integer N > 2, a polynomial f(x) = Y^k=i fk xk 
(mod N), where /i, / 2 , . . . , }k are non-negative integers and K > 1, is said to be a PP over 
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Z N when f(x) permutes {0, 1, 2, . . . , N — 1} lfT2l - |[T4ll . It is immediate that we can use this 
constant-free PP without losing generality in our quest for an inverse PP by the Lemma 2.1 
in 0. 

In this correspondence, let the set of primes be V = {2,3,5,...}. Then an integer N can be 

factored as iV = J~[ jo njVp , where p's are distinct primes. In addition, n N:P > 1, for a finite 

p&v 

number of p and n N:P = otherwise. 

Theorem 2.1 (Eft, Eft): Let N = J] p njv f and denote a divides /3 over Z by a|0. The 

pG-P 

necessary and sufficient condition for a quadratic polynomial /(x) = f\x + /2^ 2 (mod N) to 
be a PP can be divided into two cases. 

1) 2\N and A\ N (i.e., n N>2 = 1) 

/i + f% is °dd, gcd(/i, y) = 1 and / 2 = Y[ P nf ' p > n f,p — h such that p / 2 and 

P&P 

n N , P > 1- 

2) Either 2 f AT or 4|iV (i.e., ^,2 7^ 1) 

gcd(/i, N) = 1 and / 2 = n p nf > p ,rif. p > 1, Vp such that n N)P > 1. 
Theorem 2.2 ( 4701/ ): Let a, /3 be any integers and A" be a positive integer. The linear con- 
gruence ax = [3 (mod N) has at least one solution if and only if j\(3, where 7 = gcd(a, N). 
If 7 then it has 7 mutually incongruent solutions. Let x be one solution, then the set of the 
solutions is 

N 2N (j-l)N 

Xq, Xq H , Xq H , . . . , Xq H 

7 7 7 

, where x is the unique solution of ^x = ~ (mod ^). 

Definition 2.3 ( /[/j]/, $33$ ): Two polynomials /i(x) = Xlf=i /l.ft^ and / 2 (x) = J]f=i /2,fc£* 
of degree K are called congruent polynomials modulo iV if = / 2 ,fc (mod A/"), where 1 < 
k < K and equivalent polynomials modulo iV if fi(x) = / 2 (x) (mod AQ, where < x < N— 1. 

Definition 2.4 ( /O/, £75]/, /f7?l/): A polynomial z(x) = X/*~=i ^fc 2 ^ ( m °d -^0 is called a non- 
trivial zero polynomial of degree K modulo A" if zk ^ and z(x) = 0, < x < N — 1. 
Specifically, z(x) = is a trivial zero polynomial. 

Proposition 2.5 ( KT3\l . /fi?l/): If two polynomials /i(x) and ,/ 2 (x) are equivalent but not 
congruent, there exists a non-trivial null polynomial z(x) such that /i(x) — / 2 (^) = z{x) 
(mod iV). 
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Definition 2.6: Let f(x) be a PP. A PP of least degree has a least degree among all equivalent 
polynomials of f(x). 

The following proposition was proposed in [Q~3), f[i"4l . The proof is shown for its simplicity. 

Proposition 2.7 ( £T3]/, gBH): Let f(x) = J2k=i fk% k (mod N), where f K ^ and K > 
iV. Then there exists an equivalent polynomial of f(x) such that the degree of the equivalent 
polynomial is less than N. 

Proof: Let z(x) = fx • x K ~ N ■ Ylk=o( x ~~ Clearly z(x) is a zero polynomial. Let 
f(x) = f(x) — z(x), then f(x) = /(s) but deg{/(x)} < deg{/(x)}. By applying this repeatedly, 
an equivalent polynomial of degree equal to N — 1 can be found. ■ 

Proposition 2.8 ( [6]): Let f(x) = f\x + f2X 2 be a QPP and let k be an integer such that 
k > 1. Let us take /(x) such that 2 \ f t when 2|iV and 4\ N. Then / a + A;/ 2 is an unit for all 
k > 1, i.e., fi + A;f 2 is invertible and , n , is well defined. 

Proof: By Theorem I2.2L an element f\ + /c/ 2 in integer rings is an unit if and only if 
gcd(/i + kf 2 , N) = 1. We show that gcd(/i +kf 2 ,N) — 1. 

1) 2|JV and 4 { iV (i.e., n Wi2 = 1) 

In this case there exist two equivalent QPPs [9], i.e., fix + / 2 x 2 and (/i + (/ 2 + y)a: 2 , 

where 2 f /i. Let us take a polynomial /(x) = f\x + / 2 x 2 such that 2 f /i. Suppose that 
gcd(/i + kf 2 ,N) ^ 1. Then there exists a prime p such that + kf 2 ) and p|iV. By 
Theorem 12. 1[ if p\N, then p|/ 2 but p\ f\. A contradiction. 

2) Either 2\N or A\N (i.e., n Wi2 ^ 1) 

In this case, there exist one (if 2 \ N) or two (if 4|iV) equivalent QPPs [9]. In either case, 
by Theorem 12.11 and a similar argument in (1), gcd(/i + kf 2 ,N) = 1. 

■ 

Since the inverse of only one of the equivalent polynomials is sufficient for our purposes, f(x) = 
fix + f 2 x 2 such that 2\ fi will be considered in the rest of the correspondence. The following 
corollary is an extension of Proposition 12.81 

Corollary 2.9: Let fx, f 2 and iV be the integers in Theorem 12.11 and let k, kx and k 2 be 
integers such that 1 < kx < k 2 . Let us take f(x) such that 2 { fx when 2\N and 4 { N. Then 

gcd{ntfa(/i + k h),N} is an unit. 

Proof: This is a direct consequence of Proposition 12.81 ■ 



JAN. 23, 2011 



5 



III. Inverses of Quadratic Permutation Polynomials 

In this section, we derive a necessary and sufficient condition for a QPP to admit a least degree 
inverse in Theorem 13. 101 (main Theorem). We also explicitly find the inverses in Algorithm U 
This section is organized as follows. We first show that the problem of finding inverse PP(s) 
of least degree is equivalent to solve a system of linear congruences. Then we show that the 
inverses can be found by factoring the matrix for a system of linear congruences and solving it. 
We also show that solving the system of linear congruences can be much simplified and finally, 
by showing the number and the form of zero polynomials, we find all the inverses of a QPP. 

Lemma 3.1: Let f(x) = f±x + f2X 2 (mod N) be a QPP. Then there exists at least one inverse 
g(x). Further, finding all inverse PP(s) up to degree — 1 is equivalent to solving a system of 
linear congruences, 

Ag = b (mod N), 

where 

a id = (1^ + ^0,1 <i,j<N-l, 
g = [gi, 92, 9n-i\ T , and b = [b u b 2 , 6jv-i] T = [1, 2, N - if. 

Proof: Since the set of PPs forms a group under function composition, the existence of an 
inverse for a QPP is guaranteed flTl, lfl"5l . Let g(x) be an inverse PP of f(x) and suppose that 
deg{g(x)} > N. Then by Proposition I2.7[ it can be reduced to an equivalent polynomial of 
degree less than N. 

Since g(x) is an inverse, (gof)(x) = x, where < x < N — l. The equivalence of (gof)(x) = x 
and Ag = b is shown by evaluating (g o f)(x) = Y^k=i 9k(fi x + f2X 2 ) h = x at each point 
1 < x < N — 1. Note that (g o /)(0) = trivially holds. Consequently, solving Ag = b is 
equivalent to finding all the inverse PP(s) up to degree N—l. Since the number of inverse PP(s) 
up to degree iV — 1 is finite, there exists a least degree inverse. 

■ 

Lemma 3.2: Let A be an iV — 1 by N — 1 matrix in Lemma 1370 Then A = LDU, where L, 
D and UareiV — lbyiV — 1 matrices as shown below. 
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L is an N — 1 by N — 1 lower triangular matrix such that 

(;) ■ mtru + m m>j 

otherwise 

D is an TV — 1 by N — 1 diagonal matrix such that d i:i = i\, where 1 < % < N — 1. 
U is an TV — 1 by N — 1 upper triangular matrix such that 

1 if i = j 

q(*j)V(* J )r (j ) if i < j 

otherwise 



u 



■':J 



q^'i) is an 1 by j matrix such that = (ifi + i 2 1 f 2 ) k l , where 1 < k < j 
i s a j j-,y j U pp er triangular matrix such that 



y(*j) 

" nL-iW^) 

and = [0, 0, 0, 1] T is a j by 1 matrix. 
^y(fcj') i s a j by j upper triangular matrix such that 



if i = 1 
otherwise 



(M 







if m > n 



(Jfe/i + k 2 f 2 ) n - m - x otherwise 
where 1 < m,n < j . 

Proof: See Appendix A. ■ 
The factorization in Lemma I3T21 is similar to LDU decomposition except that L has not Is on 
the diagonal [[8]|. 

Lemma 3.3 ( ffiTll . KWJ): Let A, L, D and U be the matrices in Lemma 13.21 Then Ag 
b Dh = e, where h = Ug and e = L _1 b. 

Let us identify iV — 1 by 1 matrices g = [gi, g 2 , ...,gN-i] T , h = [hi, h 2 , h N _i] T with g(x) = 
J2k=i 9kX k , h(x) = J2k=i hk% k , respectively. Then the degree and the number of g and h are 
equal. 

Proof: Since all the diagonal elements of L are units by Corollary 12.91 L is an unit [UTI . 
Thus Ag = b <^ DUg = L x b. Let h be an N — 1 by 1 matrix such that h = Ug. Since U 
is also an unit, the degree and the number of g and h are equal ifTTl . lfT8l . ■ 
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In the following, two corollaries of Lemma 13.31 are shown. 

Corollary 3.4: The linear congruence Dh = e has at least one solution, i.e., there exist hk's 
such that dk,k • hk = e^, where 1 < k < N — 1. 

Proof: Suppose that for some k, there does not exist hk such that d^k • hk = e&. Then there 
does not exist a solution of Dh = e. By Lemma [3731 there does not exist a solution of Ag = b, 
which contradicts Lemma 13.11 ■ 
Corollary 3.5: Let us consider the linear congruence Ag = b. There exists a least degree 
inverse g such that deg{g} = K if and only if e K ^ and e k = 0, where K + 1 < k < N — 1. 

Proof: 
( ) 

Let g be a least degree inverse such that deg{g} = K. By Lemma [3731 the degree of h is also 
K, i.e., h K ^ and /i^ = 0, where K + 1 < < N — 1. Since Dh = e, = 0, where 
K + l<k<N — 1. Suppose that = 0, i.e., dx,^ ■ h K = 0. Let us define an iV — 1 by 1 
matrix h' such that 

h k , 1 < k < K - 1 



^1 



0, < A; < iV - 1 



Then h' also satisfies the linear congruence Dh' = e. Let g' be an iV — 1 by 1 matrix such that 
h' = Ug', then g' is also an inverse. Since deg{h'} = deg{g'} < K by Lemma 13731 g cannot 
be a polynomial of least degree. This contradicts the assumption. Consequently, e K ^ 0. 

(^) 

(1) Suppose that deg {g} > K. Then by Lemma 13731 deg {h} = deg {g} > K, where h = Ug. 
Let us define an iV — 1 by 1 matrix h' such that 

h k , 1 < k < K 

0, K + l<k<N - 1 



h' k 



Then deg{h'} = K and h' also satisfies the linear congruence Dh' = e. Then again by 



Lemma [3731 deg {g'} = K, where h' = Ug'. Since deg{g'} = K < deg {g} and Ag' = b, g 
cannot be a polynomial of least degree. 

(2) Suppose that deg {g} < K. Since g K = 0, h K = by Lemma [3731 Consequently, d KK -h K = 

eK = 0. This contradicts the assumption e K ^ 0, thus deg {g} cannot be less than K. 

By (1) and (2), deg{g} = i^. ■ 
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In Lemma 13731 since all the entries of L, Lr 1 , D, U and e can be computed for the given 
f(x) = fix + f2X 2 , finding the inverse of f(x) reduces to solving N — 1 linear congruences 
Dh = e and h = Ug. However, the cost of computation for the matrices can be substantial for 
a large N. 

The computational complexity is shown to be significantly reduced by the following lemma and 
corollary. The following lemma shows that the degree of the least degree inverse has an upper 
bound. 

Lemma 3.6 ( [7]): Let N = f\ p nN ' p . If f(x) is a QPP, then the inverse PP has degree no 

pev 

larger than max 72,7V »• 
pev ,F 

Proof: Since the set of PPs is a finite group as shown in Lemma 13 . 1 L there exists an integer 
m called an order such that the m-fold composition of f(x) with itself is an inverse PP [fl5l . 
Let f( n \x) be n-fold composition of f(x) with itself. It is shown that the coefficient of the 
degree k term of f^ n \x) is divided by ft 1 as follows. For f^(x), it is clear that / 2 divides the 
coefficient of the degree 2 term. If the coefficient of the degree k term in f^ n \x) are divisible 
by /I" 1 , then the coefficient of the degree k term in / (n+1) (x) = /i(/ (n) (a:)) + /2(/ (n) (^)) 2 are 
also divisible by ft 1 - By induction, the coefficient of the degree k term of f^ n \x) is divided 

by ft 1 - 

Suppose k > m&xn Np + 1. Since / 2 is divisible by the factors of N, Nlft 1 , i-£-, ft 1 = 
pev 

0. Consequently, there exists an inverse f^(x) that contains no terms of degree larger than 

maxnjv,,. ■ 
p&v ' y 

Corollary 3.7: Let us consider the linear congruence Dh = e in Lemma [3731 

For all k such that k > maxn Nv + 1, et = 0. 

pev 

Proof: 

Let the degree of the least degree inverse be K. Since there exists an inverse such that the 

degree of the inverse is no larger than maxnjv„ by Lemma [3761 K < maxn/vp. Consequently, 

P&V ' pev 

by Corollary 13.51 e& = 0, where K + 1 < maxn/v iP + 1 < k < N — 1. ■ 

peV 

By Corollary 13.71 only maxn Nv by maxn Nv leading submatrices (the upper-left corners of 

pev pev 

matrices) of L, L _1 , D, U and a max% iP by 1 leading submatrix of e are required to be 

peV 

computed for finding the inverse of least degree. For example, let TV = 2 18 -3 2 -5. Then maxn^ p = 

peV 

max{18,2, 1} = 18, thus only 18 by 18 leading submatrices of L, L _1 , D, U and a 18 by 1 
leading submatrix of e need to be computed instead of N — 1 by TV — 1 submatrices of L, L~\ 
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D, U and a N — 1 by 1 submatrix of e. 

In the following proposition and corollary, it is shown that the computational complexity for 
the matrices can be further reduced. 

Proposition 3.8: Let e be an N — 1 by 1 matrix in Lemma 13.31 Let also C^, where k > 0, 
be a sequence of integers known as Catalan numbers. The kth Catalan numbers are given by 

r J— ( 2k \ (2fc)! 

k k+l\k) (k + l)\-k\' 
A recurrence relation for C& is 

Cfc = — ; — ; : - • Cfe-l, k > 2, 

k + 1 

i.e., C = 1, Ci = 1, C 2 = 2,C 3 = 5, C 4 = 14, C 5 = 42, C 6 = 132 Then, 

k\ ■ C fe _! • (-f 2 ) h ~ l 



efc 



where 1 < k < 50. 

Proof: Let L, e and b be the matrices in Lemma [3731 Let also the kby k leading submatrix 
of L, k by 1 leading submatrices of e and b be L', e' and b' respectively. 
The following statement, b' = L'e', was verified to be correct for 1 < k < 50. 

k 

bk = k = ^ Ik^n ' e n 
n=l 



n=l >- v \ / m=k 



k\u 2 ::hh + mh) -(k-ny.- c ^- { h) 

Since b' = L'e', it is clear that b' = L'e'. Consequently, e' = L' _1 b' for 1 < k < 50. ■ 

Corollary 3.9: Let N = Yi p nN ' p < 2 50 and let also e be an N — 1 by 1 matrix in Lemma [3731 

pG-P 

If efe = for some k, then e n = for n > k + 1. 

Proof: Let TV = ]~J p njv >p < 2 50 , then clearly max n^ p < 50. By Lemma [3761 and Proposi- 
tionT 



efe 



n m =i (/i+ m /2) 

0, 51<A;<iV-l 
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Since /1+/2 is an unit, e\ = j-^rjq ^ 0. Suppose that e k = for some k. Since C k = 2<y2 k ^ -C k -i 
for k > 2, = 2 (2k - 1). Thus e k+1 = g^H=g^ . Consequently, e fc+1 = if e fc = 

for some k. By induction, if e k = for some k, e n = 0, for n > k + 1. ■ 
We are not aware of a closed-form expression of e when k is larger than 50. However, the 
investigation on the inverse of a QPP is not restricted under this condition since the interleaver 
size N is far less than 2 50 in practice. By Proposition 13.81 and Corollary 13.91 matrices L and 
L 1 need not to be computed for solving Dh = e. Combining Lemma I3TT1 13.21 Proposition 13.81 
and Corollary 13.91 we state the main theorem. 

Theorem 3.10 (main Theorem): Let N = Y[ p nN,p < 2 50 . The necessary and sufficient con- 
dition for a QPP to admit a least degree inverse g such that deg{g} = K is finding a smallest 
integer K > 1 such that 

(K+ 1)! ■ C K ■ = modJV 
and the number of inverse PP(s) is 

K 

IJgcd(fc!,iV). 



k=l 



Let us slightly abuse the notation in this theorem (and in examples and Algorithm U) by writing 
D, U, g, h and e for K by K leading submatrices D, U and K by 1 leading submatrices g, 
h, e, respectively. The inverse PP(s) can be found by using either (1) or (2). 

(1) Find all h's such that Dh = e and corresponding g's such that h = Ug. 

(2) Find a h such that h k = J^t'Sr/^ — 77, corresponding g and add it TTf-i gcd( k\,N) zero 

n m =i(/i+ m /2) s " 1 

polynomials. 

Zero polynomials of degree K are J2k=i | g cd(fc!,jv) ' T k ' H^o( x ~ m )\ wnere < r fe < 

gcd(A;!,A^) - 1. 

Proof: The necessary and sufficient condition is shown by combining Corollaries 13.51 and 
13.91 By Corollary 13.51 g is an inverse of least degree such that deg{g} = K if and only if 
e K and e k = for K + 1 < k < N - 1. By Corollary if e k = for some k, 
then e n = for n > k + 1. Thus g is a least degree inverse such that deg{g} = K if and 
only if eK ^ and e K+1 = 0. Since e x ^ 0, finding the degree of the least degree inverse 
is equivalent to finding the smallest K such that 6k+i = 0. Consequently, the necessary and 
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sufficient conditional for a QPP to admit a least degree inverse is (K + 1)! • Ck • f% = since 
e *+ 1 = itVM =0^(X + l)!-6^-/ 2 =0. 

The number of solutions of linear congruences Dh = e is nf=igcd(fc!,iV), since kth linear 
congruence is d kjk • h k = e k and gcd(d k iV) = gcd(&!, AT). By Lemma [331 the number of 
solutions of Ag = b is also Ylk=i g cc K&!, N). 

The complete solution set can be obtained by exhaustively solving Dh = e and h = Ug. 
An alternative is to find one solution h and g such that Dh = e, h = Ug and add it zero 



polynomials of degree K. Consider kth linear congruence Dh = e, i.e., d ktk • hk = eu- Clearly 

h k = rrS-iTfr \ \ is a solution of d k>k ■ h k = e k , i.e., k\-h k = k ^2k-i'}~{ 2)k f r - The number and 

ll m =i Ui+mfv ' Um=i ih+mh) 

form of zero polynomials are shown in Appendix B. 



IV. Examples 

We present four examples to illustrate the necessary and sufficient conditions of Theorem 13 .101 
The first and second examples consider interleavers that was investigated in @ and |fT9l . The 
third example shows the exact least degree for inverse polynomials can be less than an upper 
bound derived in and the fourth example shows the necessary and sufficient condition for a 
QPP to admit a least degree inverse g such that deg {g} = 2, 3, 4 and 5. 
All good quadratic interleavers found in Table [I] admit low degree quadratic inverses. This 
observation may not be completely surprising because [J2], [J3J shows that good interleavers 
should require the second degree coefficient to be relatively large (which works toward satisfying 
Theorem 13.101 ) but bounded by some constraints. 



1) Let f(x) = fix + f 2 x 2 mod N, where N = 1504 = 2 5 • 47, fx = 23 and f 2 = 2 • 47. The 
smallest K such that (K + 1)! • C K ■ = is 3. 

By Lemma [3T2l 3 by 3 matrices D, U and a 3 by 1 matrix e are computed as follows. 



1 










1 


117 


153 




797 





2 





,u = 





1 


539 


>e = 


188 








6 










1 




752 



'if K = l, f(x) is a linear PP. 
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Let us now exhaustively solve the equation Dh = e (mod 1504). From di t i -hi = e\, 
d-2,2 ■ h 2 = e 2 , d 3i3 • h 3 = e 3 , we can obtain hi = 797, h 2 = 94,846 and h 3 = 376, 1128, 
respectively. Since gcd(l!,iV) = 1, gcd(2!,iV) = 2 and gcd(3!, N) = 2, the number of 
solutions is 4. Let us choose hi = 797, h 2 = 94 and h 3 = 376. We obtain g 3 = h 3 = 376, 
92 = h 2 — u 23 ■ g 3 (mod 1504) = 470 and gi — hi — U\ 2 ■ g 2 — -u 13 ■ g 3 (mod 1504) = 1079 
by solving Ug = h (mod 1504). 
2) Let N = 6016 = 2 7 ■ 47, f x = 23 and / 2 = 2 ■ 47. The least degree is 4. A 4 by 4 matrix 
U and a 4 by 1 matrix e are computed as follows. 



U = 



1 


117 


1657 


1357 




3805" 





1 


539 


507 




188 








1 


1454 


,e = 


752 











1 




3008 



Let us compute h k such that h k = f£^j^ for each k. Then h = [3805, 94, 4136, 4888] r 
and g = [1831, 3854, 1880, 4888] T . 

3) Let N = 2 24 , fi = 26119 and f 2 = 2 ■ 3 ■ 41 • 179 The least degree K is 12, which shows 
the upper bound 24 obtained by the technique in is not tight^ 

4) The necessary and sufficient condition for a QPP to admit a least degree inverse g such 
that deg{g} = K = 2, 3, 4, 5 is 12/f = 0, 120/f = 0, 1680/ 2 4 = and 30240/f = 0, 
respectively. This formula is also shown in 01, [O for K = 2 and in 01 for K = 3. 



188 QPP based interleavers have been proposed in 3GPP LTE [19]. Most of the interleavers 
proposed in [19 ] admit a quadratic inverse with the exception of 35 of them. In Table [III all of 
the interleavers that do not admit quadratic inverses are listed with their respective inverse PPs 
of least degree computed using Algorithm [0 

2 An inverse g(x) is 7612343a;+4897586a; 2 +352440a; 3 +2867432a; 4 + 13756448a; 5 + 13890368a; 6 +915200a; 7 +2679424a; s + 
6846976x 9 + 5217280x 10 + 53248a; 11 + 1478656a; 12 . 
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TABLE I 
Algorithm 1 



An algorithm for finding the inverse PP(s) of least degree for a QPP f(x) = f\x + /2a; 2 (mod N) 

1. If 2\N, 4 \ N and 2|/i, let f(x) be such that f(x) = (/1 + f)x + (/ x + f )a; 2 . 

2. Find the smallest integer K > 1 such that (A' + 1)! • • = 0, where C = 1 and C fc = ^ ( 2 fc fe ). 
Then, the least degree of the inverse PP(s) is K. 

3. Compute K by K matrices D, U in Lemma [3~2l and K by 1 matrix e in Proposition 13.81 

4. There exist two methods for finding the solution set of Ag = b Dh = e, h = Ug. 

(1) All the h's and g's can be found by solving K linear congruences Dh = e and h = Ug. 

g's can be computed by by back-substitution. 

Note that g K = h K and g k = h k - J2m=k+i ' 9m for 1 < k < K - 1. 

(2) Find one inverse and add it zero polynomials of degree K. 

Compute hk = — T) ^ — ^ — ^ anc ' corresponding g such that h = Ug. 
Convert K by 1 matrix g into a polynomial and add it z(x) = Ef=i {^TjZ. 
where 1 < r fe < gcd(fc!, N) - 1. 



V. Conclusion 

We derived in Theorem 13.101 a necessary and sufficient condition to determine the least degree 
inverse for a QPP. We also provided an algorithm to explicitly compute the inverse PP(s). 
188 QPP interleavers were proposed in 3 GPP LTE |fl9l . Most of the QPP interleavers in [[191 
admit a QPP inverse. We applied the theory in this correspondence to tabulate all inverse PPs 
of degree larger than two. Further, it was shown that inverses of good interleavers in [[191 have 
low degrees and a possible explanation is given. 
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TABLE II 

Inverse PPs of Least Degree for 3GPP LTE Interleavers without Quadratic Inverses 



length 


OPP 


An Inverse PP of Least Degree 


length 


OPP 


An Inverse PP of Least Degree 


928 


15x + 58a; 2 


31a; + 290a; 2 + 232a; 3 


4544 


357a; + 142a; 2 


4509a; + 994a; 2 + 2840a; 3 


1056 


17a; + 66a; 2 


1025a; + 726a; 2 + 792a; 3 


4672 


37a; + 146a; 2 


2557a; + 1022a; 2 + 4088a; 3 


1184 


19a; + 74a; 2 


779a; _|_ 74a; 2 _|_ 296a; 3 


4736 


71a; + 444a; 2 


2935a; + 3996a; 2 + 3552a; 3 


1248 


19a; + 78a; 2 


427a; + 78a; 2 + 936a; 3 


4928 


39a; _|_ 462a; 2 


1927a; + 1078a; 2 + 616a; 3 


1312 


21a; + 82a; 2 


781a; + 574a; 2 + 984a; 3 


4992 


127x + 234a; 2 


511a; + 2730a; 2 + 4056a; 3 + 2184a; 4 


1376 


21a; + 86a; 2 


557a; + 602a; 2 + 344a; 3 


5056 


39a; + 158a; 2 


3079a; + 2054a; 2 + 1896a; 3 


1504 


49a; _|_ 846a; 2 


353a; + 282a; 2 + 376a; 3 


5184 


31a; + 96a; 2 


3679a; + 1632a; 2 + 1152a; 3 


1632 


25a; + 102a; 2 


1273a; + 306a; 2 + 408a; 3 


5248 


113a; + 902a; 2 


2833a; + 410a; 2 + 4264a; 3 + 328a; 4 


1696 


55a; + 954a; 2 


663a; + 530a: 2 + 424a; 3 


5312 


41a; + 166a; 2 


3401a; + 498a; 2 + 3320a; 3 


1760 


27a; + 110a; 2 


163a; + 990a; 2 + 1320a; 3 


5440 


43a; + 170a; 2 


1107a; + 1530a; 2 + 680a; 3 


1824 


29a; + 114a; 2 


1541a; + 1710a; 2 + 1368a; 3 


5504 


21a; + 86a; 2 


2621a; + 5074a; 2 + 1032a; 3 + 5160a; 4 


1888 


45a; + 354a; 2 


21a; + 1534a; 2 + 472a; 3 


5568 


43a; + 174a; 2 


1651a; + 1566a; 2 + 4872a; 3 


1952 


59a; + 610a; 2 


579a- + 1586a; 2 + 488a; 3 


5696 


45a; + 178a; 2 


3829a; + 5518a; 2 + 4984a; 3 


2112 


17a; + 66a; 2 


1025a; + 1782a; 2 + 792a; 3 


5824 


89a; + 182a; 2 


409a; + 3458a; 2 + 3640a; 3 


2944 


45a; + 92a; 2 


1701a; + 1748a; 2 + 2208a; 3 


5952 


47a; + 186a; 2 


95a; + 930a; 2 + 3720a; 3 


4160 


33a; + 130a; 2 


3057a; + 1430a; 2 + 1560a; 3 


6016 


23a; + 94a; 2 


1831a; + 3854a; 2 + 1880a; 3 + 4888a; 4 


4288 


33a; + 134a; 2 


3281a; + 1474a; 2 + 2680a; 3 


6080 


47a; + 190a; 2 


2943a; + 950a; 2 + 2280a; 3 


4416 


35a; + 138a; 2 


347a; + 2346a; 2 + 552a; 3 









Appendix 

(A) [Lemma Q 

We use two-fold induction and prove A = LDU by showing that column-reduced form of A 
is equivalent to LD. 

Let us define an A — 1 by A — 1 elementary matrix such that 

1 if m = n 



t 



(m) 



-Ui 



Hj if m = i, n = j , 

otherwise 
where 1 < % < j - 1, 2 < j < A - 1 and 1 < m, n < A - 1. 

Let T = T (1,2 ) • T^ 1 - 3 ) ■ T( 2 ' 3 ) ■ ■ ■ T^ N ^ ■ ■ ■ T^- 2 -^ 1 ), then it is easily verified that T = IT 1 . 
Let us also define N — 1 by N — 1 lower triangular matrices L^'^ such that 

L(ij) _ at (1 ' 2) ■ T (1,3) • T (2,3) ■ ■ ■ T (1,i) • ■ ■ T (i ~ 1,J ' ) T (ij ' ) i.e., 
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LO'-ijOt^J) if i = 1 
L(*-y)TW) if 2 < i < j - 1 



Since U is an unit, A = LDU if and only if 

att -1 = AT = AT (1 ' 2) • T (1>3) • T (2 ' 3) • • • T (i,j) • ■ ■ T (1,JV-1) • ■ ■ T^ -2 '^ - ^ 
s v ' 

= L^" 2 ^- 1 ) = LD. (1) 
We use induction on j and prove eq. © by showing that L^' -1 '^ is as follows. 

{n! • l m>n if 1 < n < j 
(2) 
a m ,n if j + 1 < n < N - 1 

Upon completion of column reduction, j = N — 1, thus eq. (OQ) holds. 
We first show that eq. © holds for j = 2. 

By definition, L (1,2 ) = AT^ 2 ). Since tj* 2 2) = -u 1(2 and w 1>2 = q^v^M 2 ) = fc + f 2 , 

^,2 — ~ u l,2 " Sm,l + a m,2 

= -(/i + / 2 )-(m/ 1 +m 2 / 2 ) + (m/ 1 + m 2 / 2 ) 2 

= (m/i + m 2 f 2 ) ■ {{m - l)h + (m - l)(m + l)/ 2 } 

= ( m -l). m .(/ 1 + / 2 ).{/ 1 + ( m + l)/ 2 } 



Consequently, 



,(1,2) 
l m,2 



if m — 1 



2! • / m>2 if m > 2 

Thus eq. © holds for j = 2. Suppose now that eq. © holds for j > 2. For each j, we use 
induction on % and show that eq. © holds. Upon completion of induction on i, we show eq. © 
holds for j + 1. 

i 

S+l = [Il( m - Wl + ( m + ■ q^^V^^M^ 1 ) (3) 

fc=0 
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In the following, h^ +1 > is shown below in matrix form. 



low\ 
low 2 

loiVi 

low i+1 



loWj 



Iown- 



coli 

( v..i 1A 

1! • f 2)1 

1! ■ k,i 
1! • k+x,x 



col 2 


2! • h,2 

2! ■ k,2 
2! • k+1,2 



2\-l 



J-2 



1! " ki 

\ 1! ■ In-1,1 2! ■ In-1,2 



colj 


C0lj + i 














(X\ JV— 1 








a 2j+2 


02,JV-1 






'ij'+l ~ U 
lW+1) 


a i,j+2 
a i+l,j+2 


Oi+l,iV- 






a jd+2 


a j>N -i 



ajv-i,jv-i y 



The elementary matrix T ( - 1 ' 3 ' +1 - ) subtracts w x j +1 times column 1 from column j + 1 of L^ -1 '^. 
We show that L^ -1 '^ multiplied by T^ 1 leaves other columns unchanged except the column 
j + 1 and creates a zero in the (1, j + 1) position of = L^ _1 ^T( 1 « 3 " +1 ). 

When z = 1, eq. © holds, since 



I 



(iJ+i) 

mj'+l 



-«l,i+l • I' • ^m,l 

_ q (U+i) V^+Dr^ 1 ) ■ W (/1 + m/ 2 ) + (m* + m 2 / 2 ) J+1 
-(/1 + / 2 ) J • {mh + m 2 / 2 ) + (mh + m 2 f 2 y +1 
(mh + m 2 f 2 ) ■ Umh + m 2 - (f ± + 



mfi + m 2 f 2 ) • |(m - 1)/! + (m - l)(m + l)/ 2 j • ^ |(m/i + m J / 2 ) 



(/1 + 



fc=0 



m(m 



- !)(/! + m/ 2 ){/! + (m + l)/ 2 } • c£ m > j+ VW^ +1 W + V 



[U(m-k){f 1 + (m + k)f 2 } 

k=0 



q 



(m,j+l)y(2,i+l) r (3+l) i 



Thus l<tf£> 



as desired. 

Suppose now that eq. © holds for i. T^ +lj+1 ^ subtracts iij+ij+i times column i+1 from column 
j + 1 of L^ ,J-+1 s where 2 < z < j. In the following, it is shown that that L/ M+1 ) multiplied by 
leaves other columns unchanged except the column j + 1 and creates a zero in the 
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(i + + 1) position of = L^'+^T^ 1 '^ 1 ). 



m,j+l 



' L m,i+1 L m,j+1 

• (i + 1)! • / m ,i+l + *m5+l 



ra+j 
fc=m 

i 

[J|(m - k){fi + (m + fc)/ 2 }] • q (^'+ 1 )v (i+lj+1) r ( - ? ' +1) 

fc=0 

j 

- _ q (<+iJ+i)v( i + 1 J'+ 1 )r^ +1 ) • [fj(m - + (m + fc)/ 2 }] + 

fc=0 

i 

[J|(m - fc){/i + (m + fc)/ 2 }] • q(^+i)v( i+1 ^ +1 )r^' +1 ) 

k=0 

i 

= [ J|(m - fc){/i + (m + fc)/ 2 }] • {q(^ +1 ) - q (<+iJ+i)}v( i+1 ^ +1 )r^ H 



i) 



q 



k=0 



= 1, m/x + m 2 / 2 , (Wi + m 2 f 2 ) j \ [l,(i+ l)f\ + (i + l) 2 f 2 , {(i + l)h + (i + 1) 2 / 2 } J ] 
= {m + l)}{/i + (m + i + l)/ 2 } • 

0, 1, m(/i + m/ 2 ) + (< + l){/i + (i + l)/ 2 }, + ^/2)} J_1_n {(^ + + (i + l )h)Y 

n=0 

= {m - (i + 1)}{/! + (m + i + l)/ 2 } • q (^+i)w (i+lj ' +1) . 
Thus, 



[f](m - + (m + k)f 2 }] ■{m-{i+ \)}{f\ + (m + i + l)/ 2 } 



fc=0 



q(mj+l) w (i+lj+l)y(i+lj+l) r (j+l) 



i+1 



= Y[(m - k){h + (m + k)f 2 } ■ q (^ +1 )v (i+2j+1) r (i+1) 



fc=0 
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Consequently, ^^j+i — and eq. © holds for i + 1. 

We now show that eq. © holds for j + 1. Let % = j in eq. ©. Then 

q (™J+l)y(j+l,j+l) 

0, 1, ...] • W (i ~ lj ' +1) W (j ~ 2 ' j+1) • . . w (1,j+1) 

o,o,i,...] . w^- 2j ' +1 )---w( lj ' +1 ) 



. w (iJ+1) W ( ^ 1J+1) W (i ~ 2,j+1) ■ . . w (1,j+1) 



[0,0,...,!]. 



Thus q( m J+ 1 )v^' +lj ' +1 M' ? ' +1 ) = 1. Consequently, 1%*$ = 0, where m < j and 



fc=0 



(i + i)! 



/77 



^ k=m 



= {j + 1)\ ■ l m>j+U 

where j + l<m<N — 1. Consequently eq. © holds for j + 1. 
(B) [The number and the form of zero polynomials of degree K] 

We show the number and the explicit form of zero polynomials of degree K, where K < N — 1[] 
In Lemma lATTl the necessary and sufficient conditions for a polynomial to be a zero polynomial 
is shown and in Lemma IA.2I and IA.31 the number and the explicit form of zero polynomials of 
degree K are derived by using Lemma IA.1I 
Let us define z n (x), where < x < jV — las follows. 



z n (x 



z n -\(x + 1 mod N) — z n _i(x mod N) 



Lemma A.l: The following statements are equivalent. 
(1) z{x) = 0, where < x < N - 1. 



if n = 

if 1 < n < K 



3 A different proof is shown in 1161 pp. 245] and [13] for the explicit form of zero polynomials. 
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(2) z(x) = 0, where < x < K. 

(3) z n (0) = 0, where < n < K. 

Proof: 
( (1) (2) ) 
Trivial. 

( (2) (3) ) 

It is easily shown by induction that if z n (x) = 0, where < x < K — n, then z n+ i(0) = 0, 
where < x < K — (n + 1). Consequently, (3) holds. 
( (3) =s- (1) ) 

Suppose that z n (0) = 0, where < n < K. Since zk{x) is a constant, if zk(0) = 0, then 
z K [x) = 0, where 1 < x < N — 1. 

Consider zjj-(x) = + 1) — Since z K ^i(0) = and 2^(0) = by assumption, 

<2ft:-i(l) = zk(0) + -2a'-i(0) = 0. Then by induction on x, it is shown that zk~i(x) = for 
1 < x < N - 1. 

The induction outlined above are then repeated for n = K — 2,K — 3, 2, 1. Hence, (1) holds 
as desired. ■ 
Lemma A.2: The number of zero polynomials of degree K is OfcLi gcd(A;! , N). 
Proof: 

Let A be a K by if leading submatrix of A in Lemma [3T21 and let / l5 f 2 be 1 and respectively. 
Let also L, D, tj, h and z be the corresponding leading submatrices of L, D, U, h and a zero 
polynomial of degree K. Then z(x) = 0, where < x < K is equivalent to Az = 0, where 
dij = V (mod N), z = [zi, z 2 , zk] t = U _1 h and is a by 1 zero matrix. This is shown 
by evaluating z (x) =0 at each point 1 < x < N — 1. 

In Lemma 13 .31 it is shown that L and U are units. It also holds for L and tj since all the 
elements of L and U on the diagonal are Is. Then by Theorem 12.21 and Lemma I3~3l the number 
of zero polynomials of degree K is the number of solutions of Dh = 0, i.e., rifcLiS cc K^!> ^0- 
The set of solutions of Ag = b Ag + z = b), where g has degree K, is therefore composed 
of one particular solution g and zero polynomials of degree K. 
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Lemma A. 3: Zero polynomials of degree K are of the form 

k „ k-l 

z ( x ) = £ { drA ., ^ • ^ • n ( x ~ m )}' where ° - ^ - gcd iv ) - l 

fc=l ° ^ ' m=0 

Proof: 

Consider Kth linear congruence of Dh = in Lemma [AT2l i.e., dx,K * ^jc = 0. Since = K!, 
by Theorem 12.21 and Lemma 13 .31 hx = zk = gcd ^\ n) ' Tk > wnere < tr- < gcd(-K1, A/") — 1. 
Suppose now that z^ K \x) is a zero polynomial of degree K, then z^(cc) = gcd ^, N ^ -r K - z'(x), 
where z'{x) is a monic polynomial of degree K. 

Let z'(x) = nm=o( x - m ) and consider z^ K \x) = gcd{ K\,N) ' r ^ ' Hm=o( x -m).It is clear that 
z (K)( x) = 0> where o < x < Further = -^-^-r^K! = jV-r x - i ^-^ = 0, 

thus z^ (x) = 0, where < x < K. Consequently, by eq. (2) in Lemma IA.ll (x) is a 
zero polynomial of degree K. Since r m ^ r n , where m ^ n, z^ K '(x)'s are equivalent but not 
congruent polynomials. 

Let us now consider (K — l)th and Kth linear congruences of Dh = 0, where Jik-i,k-i = 
z K -i = ged (( K / ^ 1 ) ! jv) ' n K-i and Jik,k = 0. By using a similar argument above, it is shown that 
z^ K ~ l \x) = gcd (( J ^ 1 ) !|iV ) • tk-i ■ ]lm=o( x ~ m ) is a zero polynomial of degree AT - 1. Note that 
^^(x) + z^ K ~ l \x) is also a zero polynomial of degree K. 

Applying this repeatedly for k = K — 2, K — 3,..., 2, the desired result follows, i.e., zero 

polynomials of degree K are of the form J2k=i z(k) ( x ) = Ef=i { gc d(l\,N) " T * • it^o^ ~ m ) } C 

■ 
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